Zero Temperature Hysteresis in Random Field Ising Model on Bethe Lattices: 
approach to mean field behavior with increasing coordination number z. 
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We consider the analytic solution of the zero temperature hysteresis in the random field Ising 
model on a Bethe lattice of coordination number and study how it approaches the mean field 
solution in the limit z ^ oo. New analytical results concerning the energy of the system along the 
hysteresis loop and first order reversal curves (FORC diagrams) are also presented. 
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The Random Field Ising model with zero temperature 
(T = 0) metastable dynamics is a prototype lattice model 
for understanding the dynamics of disorder driven first- 
order phase transitions and hysteresis in several complex 
systems 1]. The model was introduced more than 10 
years ago and studied analytically in the mean-field limit 
and numerically on finite dimensional lattices (3d, 4d, 
etc). It predicts the existence of a critical point in the 
systems's response to a slowly and smoothly varying ap- 
plied field. Theoretical interest in the model increased 
after many of its properties were obtained analytically 
on Bethe lattices 2] . It was found that the critical point 
exists only on lattices of coordination number z > A. In 
recent years the analysis on Bethe lattices has been ex- 
tended to obtain additional results like the trajectories 
of the first-order reversal curves in the field vs. magneti- 
zation [H — m) diagram 3] , the behavior of the different 
energy terms in the hamiltonian 4] and the demagne- 
tized states 5J. In this Brief Report we present new 
developments along three lines: (i) relationship between 
the mean field solution of the model and its solution on 
Bethe lattices of large coordination numbers, (ii) a dif- 
ferent argument for the computation of the energy terms 
that allows us to rewrite the existing results for a z = 4 
lattice in a more compact and transparent form valid for 
an arbitrary value of and (iii) the computation of the 
FORC-diagrams that allow a compact description of the 
properties of the First-Order Reversal Curves (FORC). 

The RFIM model in the mean field limit of infinetely 
weak but infinitely long range pair interactions is de- 
scribed by the Hamiltonian, 
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Here J is a ferromagnetic exchange interaction of or- 
der unity, {Si = ±l;i = 1,2,...A^} denote Ising spins, 
{hi} are independent, identically distributed, on-site 
quenched random fields, and H is an externally applied 
uniform field. The quenched fields are drawn from a 
Gaussian distribution (j){hi) of variance unity and mean 



value zero. The sums over i run over all sites of the sys- 
tem, and the sum over j runs over the entire — 1 spins 
of the system that interact with the spin Si. The factor 
N — 1 dividing J ensures that the energy of the system 
is an extensive quantity. 

The model on a Bethe lattice (deep interior of an in- 
finite Caley tree of coordination number z) is based on 
short range interactions. It is characterized by the Hamil- 
tonian, 
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J2s.S,-aJ2hA-Hj2Si (2) 



The difference between the two look alike Hamiltonians 
lies in the sum over j. In Eq. [3 the sum over j runs 
over z nearest neighbors of each site i {z = 2d^ for a 
d-dimensional cubic lattice). The factor z dividing J in 
Eq. 121 ensures that the energy remains extensive in the 
limit z ^ N. Although equations Eqs. ^ and [21 appear 
to have an identical form if z = N — 1^ but the topology 
is different in the two cases. In mean field theory every 
spin interacts with every other spin but this is not the 
case on a Bethe lattice. 

In both cases (MF and B) the local force on a site can 
be writen in terms of a generic local magnetization rrii 



Fi = Jrrii + ahi 



H 
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where in the MF case mf ^ = m^^ = ^^^^ Sj/N (ex- 
cept for a negligible small correction in the thermody- 
namic limit) and for the B case mf = YTj^i^j / ^' To 
fully specify the metastable T = dynamics one must 
fix the initial state (for instance. Si = —1 d^i H = —oo) 
and adiabatically sweep the field relaxing the spins 
according to the rule Si = sign{Fi). In the MF descrip- 
tion, since rui is independent of i, the evolution of the 
magnetization m^^{H) is trivially determined by the 
two coupled equations. 
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where the first equation relates the magnetization with 
the probability of finding a spin up and the second 
one expresses that probability as an integral of the distri- 
bution of random fields . The above equations^and[51ad- 

mit one fixed point solution for all fields if a > ac = \J~^J 
and three fixed point solutions (in a certain field range) 
below (Jc. Two of the solutions are stable and one is 
unstable. The two stable solutions are obtained by nu- 
merical iteration starting from the initial state m = — 1 
and m = 1 respectively and correspond to the two halves 
of the hysteresis loop in increasing and decreasing ap- 
plied field. Note that there is no hysteresis if <j > (jc. 
The critical point {He = 0, dc) corresponds to a non- 
equilibrium critical point of the system ^ The Bethe 
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FIG. 1: Phase diagram of the RFIM on bethe lattices as a 
function of increasing z, compared with the MF behaviour. 
The inset shows the details around the z=4 critical point, 
(color online) 

lattice allows short range fiuctuations in the environment 
of each site. Therefore rrii is not homogeneous and can 
take z-\-l different values (—1, (— z + 2)/z, • • •, {z — 2)/z, 
1). Let us introduce a variable k = 0^ - - ■ z indexing the 
environments so that = {—z + 2k)/ z. To find the 
behaviour of the magnetization (for instance, along 
the i7-upwards branch of the hysteresis loop) one should 
treat separatelly each of these environments. The cou- 
pled equations become slightly more involved 2]: 



P 



2p^ 



k=0 



1 



(6) 

+l|fe) (7) 



where P{Si = +l|/c) is the probability that Si = 1 given 
that k of its neigbours are up, which is given by 



P{Si = ^l\k) 



J{2^-l)-H 
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FIG. 2: Magnetization versus field behaviour for a j J — 0.6 
and increasing values of z as indicated. The lines are com- 
pared with the MF behaviour, (color online) 



and P* is the solution when n 
recurrence relation: 



oo (fixed point) of the 
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This means that P* satisfies: 



k=0 



k 



[PT[^-PT~ P{Si = ^i\k) (10) 



The physical meaning of P* is the probability that, along 
the up-field branch of the hysteresis loop, a spin is +1 
given that a neighbour is forced to be down. For the 
values of a and H for which it displays multiple fixed 
points only the stable fixed point obtained starting from 
P^ = will have physical meaning. For small values of 
z the results of the Bethe lattice are in striking contrast 
with those of the MF case due to the neglect of envi- 
ronment fiuctuations in this last case. The critical point 
(Hc^cTc) is absent on lattices with z = 2 (Id model), as 
well as z = 3. For z > 4, there is a value of the applied 
field H < J/z where the magnetization jumps discon- 
tinuously if <j < (jc. The size of the jump reduces with 
increasing a and vanishes as a approaches (Jc, and H ap- 
proaches He = J/z. For cr > (Jc, there is hysteresis but 
no discontinuity in the magnetization. The critical point 
for z = 4 is located at He = J/z and ac ~ .445315 J. Fig. 
shows the phase diagram for different values of z. The 
critical points for z = 4, 5 ... 35 are indicated by filled 
circles. The discontinuity in the magnetization occurs on 
a field Hdis (indicated by the thick lines) and coercivity 
(m = 0) occurs on Hcoe (thin continuous lines). The in- 
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set shows a detail of the z = 4 case, reveahng that below 

Although, the most important differences between the 
MF and B approaches are seen for small values of it 
is of interest to study the Bethe lattice hysteresis loops 
for increasing z. In the limit z ^ oo, we may expect 
the Bethe lattice results to approach those of the MF 
case. The reason is that in this limit the fluctuations in 
the exchange field at a site approach zero, and therefore 
the model on a Bethe lattice approaches the mean field 
model. This argument is not entirely transparent because 
unlike the mean field model, the nearest neighbors of 
a site on a Bethe lattice are not nearest neighbors of 
each other for any z. However, as we see from Fig^ the 
critical point in the limit z ^ oo tends to the critical 
point of the MF theory he = 0,crc = Y^2/7r. Indeed, 
Fig|21 shows that the entire magnetization curve tends to 
the MF result in the limit z ^ oo. Thus the expectation 
that the B results fall over the MF results in this limit is 
indeed born out by comparing the two numerically. The 
equivalence can also be shown analytically from Eq. EJas 
follows: In the limit z ^ oo, 

dkSiP* - -) m)dhi 



Let us now focus our attention of the exchange inter- 
action along the hysteresis loop on the Bethe lattice. In 
a previous calculation Q the exchange interaction was 
computed from the two site probability p{Si^Sj). Here 
we will show that it can be written as a single site equa- 
tion. Let us start from Eq. (11) in Ref. that gives the 
correlation of two neigbouring sites: (SiSj). By intro- 
ducing the results of Eq. (16), (21), (22) and (23) of the 
same reference one can write: 



{SiSj) = 1 -4P* +4P*Q* 



(12) 



where Q* is given by: 



fc=0 ^ 
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P{Si = +l|fc + l) 
(13) 

Similarly to P*, Q* can be read as the probability that 
(along the increasing field branch) a spin is +1 given 
that a neigbour is forced to be up. Equation (|T^ is still a 
two site equation since its last term in (|T^ contains the 
double sum with two indexes reflecting the state of the 
z—1 spins (different from Si and Sj) in the neigbour hoods 
of Si and Sj. Now we make use of the identity: 
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Taking the derivative with respect to P* one gets 



But (2P* - 1) ^ (2j9^ - 1) = in the limit z ^ oo. 
Thus the equation determining the magnetization in an 
applied field H becomes the same in the two strategies. 
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From Eq. rearranging the indices (j = /c + 1) it is 
easy to obtain: 

P*Q* = Q [P't [1 - ^*]'"^ = +11^^) 

fe=0 ^ ^ ^ 

(16) 

Introducing p5|) and IjlBfl in Hl'2|) . and taking into ac- 
count that {Ue)/N = -j{SiSj) one gets: 



{Ue)/N=--J + 



(17) 



FIG. 3: Behaviour of the reduced hamiltonian /JN as 
a function of H/ J for a / J — 0.6 and increasing values of z, 
compared with the MF behaviour, (color online) 



+ E (fc) - pr-''^-^ [1 - p{Si = +i\k)] 

k=0 ^ ^ 

This equation shows that the exchange energy can be 
computed as the energy corresponding to the saturated 
state plus an exces energy (2J/z) associated to each bro- 
ken bond. This computation also allows to write the 
full average hamiltonian as a constant term plus a 
sum over the state of a single site environment. Figure IHl 
shows the behaviour of 71/ J N as a function of the exter- 
nal field H/ J for cr/J = 0.6 and increasing values of z as 
indicated. As can be seen the behaviour also tends to the 
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FIG. 4: First-order reversal curves and FORC diagram for z=4 and three different values of cr/J, as indicated. Arrows indicate 
the value of the coercive field Hcoe/ J on the two axes, (color online) 



MF behaviour shown by a dashed line. It is interesting 
to note that even the unphysical states in the MF curve 
are recovered from the limit of the unphysical states of 
the Bethe lattice with finite z. 

The third analysis that we want to present is that 
of the FORC-diagrams. Such diagrams were introduced 
1^ in order to simplify the description of the collection 
of first-order reversal curves which describe the magne- 
tization m{H2,Hi) obtained starting from saturation, 
adiabatically decreasing the field until Hi and subse- 
quently increasing the field up to H2. The FORC- 
diagrams are computed by evaluating p = d^m/dHidH2. 
This second derivative is represented as a function of 
Hu = {H2 + Hi)/2 and = {H2 - Hi)/2. The fact 
that FORC can be computed exactly on Bethe lattices 
, which for z > 4 exhibit a disorder induced phase tran- 
sition, allows a better understanding of some interesting 
features of the FORC-diagrams. 

Fig. ^exhibits the FORC-diagrams corresponding to 
z = A. They have been numerically computed by evalu- 
ating the function p in steps of Ai^i = Ai^2 = 0.01 J/z 
for three values of a as indicated. The first interesting 
feature, which is clearly seen in Fig. ^c) is the existence 
of a well defined peak for Hu — and He — Hcoe — J/ z. 
This indicates that the maximum variation in the slope 
dm/dH2 occurs around Hi = —Hcoe and H2 = Hcoe- 
The second remarkable property is the fact that p = 
almost everywhere for He > J/z. (The only exception 
is the infinitessimally thin "ridge" entering such a region 
for (J < (Jc). The reason is that on the Bethe lattice, 
FORC become independent of H2 when H2 — Hi > 2 J/z 
since, as pointed in Ref. 3], the solutions joint (and 
merge) the main hysteresis loop. The third feature to 
notice in Fig. ^is the "valley" along an approximate line 



He — Hu ~ J/z. The reason for this valley is that the 
slope of the FORC dm/dH2 increases with decreasing 
JYi, is maximum when the reversing field Hi is close to 
—Hcoei and decreases again for more negative reversing 
fields {Hi < —Hcoe)' This valley is smooth for large a 
but becomes sharper when a ^ dc. Below ac it trans- 
forms into a discontinuity "cliff" , due to the occurence of 
the discontinuity in the hysteresis loop. The discontinu- 
ity in the loop also explains the existence of the infinites- 
simally thin ridge, when a < dc and J/z < He < Hdis- 
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